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We investigate a solvable model for energy conserving non-equilibrium steady states. The time- 
reversal asymmetry of the dynamics leads to the violation of detailed balance and to ergodicity 
breaking, as manifested by the presence of dynamically inaccessible states. Two such systems in 
contact do not reach the same effective temperature if standard definitions are used. However, we 
identify the effective temperature that controls energy flow. Although this operational temperature 
does reach a common value upon contact, the total entropy of the joint system can decrease. 
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What happens when systems in non- equilibrium steady 
states are brought into contact? For contact between 
equilibrium systems, applying the second law of thermo- 
dynamics requires them to be closed or isolated, with a 
fixed total energy. When energy is allowed to flow be- 
tween systems that are otherwise isolated from the en- 
vironment, entropy maximization determines the final 
equilibrium state and leads to temperature equalization. 
Here, we study what happens when systems in far from 
equilibrium steady states are connected and allowed to 
reach a mutual non-equilibrium steady state. 

Typically, non-equilibrium steady states are realized 
in driven systems, receiving a perpetual supply of energy 
which is removed by dissipation or contact to a thermal 
bath. These are clearly not isolated, making interpret- 
ing questions about contact between them more difficult. 
Nonetheless, one may question to what extent may non- 
equilibrium steady states be described by effective tem- 
peratures determining their properties [if, Q . 

The question of contact between far from equilibrium 
systems is motivated by striking differences between equi- 
librium systems and open, dissipative systems in contact, 
such as the breakdown of energy equipartition in driven 
granular gas mixtures [3|]. Additionally, the direction of 
energy flow and possible equalization of effective tem- 
peratures are of particular relevance in attaching a ther- 
mometer to a non-equilibrium system [4|, [5|, |6|, Lfl] . 

Soft matter systems comprised of macroscopic grains, 
colloids, foam bubbles, or living cells undergo irreversible 
dynamics. Contact between such driven dissipative sys- 
tems differs from the equilibrium situation in two aspects: 
(i) the microscopic dynamics are irreversible, thus violat- 
ing detailed balance (DB), and (ii) energy flows into the 
systems due to driving and out by dissipation [8[. As a 
step toward obtaining a theoretical understanding of this 
complicated situation, we focus on the contact between 
isolated non-equilibrium systems undergoing microscopi- 
cally irreversible but energy conserving dynamics [9] . We 
furthermore expect that this concept of an isolated non- 



equilibrium system may prove useful when studying a 
slowly degrading open system on short time scales. 

In the first part of this Rapid Communication, we in- 
troduce and study a stochastic model for isolated sys- 
tems which achieves a non-equilibrium steady state. The 
model exhibits a dynamics-induced hole of inaccessible 
states. This results from DB violation, which we quan- 
tify by means of the total probability current. The sys- 
tem is furthermore characterized by moments of the en- 
ergy distribution, which we relate to various definitions 
of effective temperatures. 

In the second part of this Rapid Communication, we 
consider contact between two such systems A and B. For 
weak coupling between them, they affect each other only 
by partitioning energy between them, with the total en- 
ergy Ea + Eb being strictly fixed. In this limit, our 
model satisfies the zeroth law of thermodynamics: if two 
systems are in mutual steady states with a third system, 
they will also be in a mutual steady state one with the 
other. This is not the case for open driven dissipative 
systems. 

Effective temperatures generally differ between the 
connected systems (as in [2, 3]). We are able, however, to 
identify an effective operational temperature 0, 0] that 
determines energy flow and equalizes at contact. This 
operational temperature is a property of an individual 
system, and does not depend on the coupling details be- 
tween the systems. This notwithstanding, in contrast to 
the second law of thermodynamics, entropy defined in 
the standard way is not maximized upon contact. 

Our model consists of N particles with positive ener- 
gies \ej\i<i<N, and is a modification of the open model 
of PJJ, L!Jj • At each time step three particles i ^ j ^ k 
in the system are chosen at random. The first two (i, j) 
exchange energy "dissipatively" with a constant restitu- 
tion coefficient < a < 1, and the third (k) receives the 
energy dissipated from them. The non-dissipated energy 
a(ei + ej) is randomly redistributed between the first two. 
Thus, the particle energies after the interaction are: 
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FIG. 1: (Color online) Configuration space occupation for 
N = 3, a = 0.6. Colorbar indicates the log of the occupation 
probability density. The triangular hole in the middle corre- 
sponds to inaccessible states. Inset: Distance to equal energy 
state. Dashed line is the theoretical prediction of Eq. (|2]). 



where < z < 1 is drawn randomly at every interaction 
from a uniform distribution, and the energies of all other 
particles in the system remain unchanged. We note that 
when a = 1, the system is in equilibrium. For a<l, al- 
though the total energy in an isolated system E = Ytf =1 ei 
remains exactly constant, its steady state is far from equi- 
librium because of the time-reversal asymmetry of the 
dynamics. We shall demonstrate this through the sys- 
tem's non-ergodicity and DB violation [l2| ]. 

For a < | not all energetically allowed states are dy- 
namically accessible. If the system starts in an inac- 
cessible state, the dynamics will never return it there. 
Nonetheless, the system has a unique steady state and all 
accessible states are dynamically connected: the system 
can evolve with time between any two of them, and en- 
semble averages coincide with time averages. For simplic- 
ity of visualization, we first consider a system of N = 3 
particles. For a = 1 our model is equivalent to an equi- 
librium system with a flat single particle density of states 
(DOS) [13], and we observe a uniform occupation of the 
constant energy surface in configuration space. For a < 1 
we find numerically a non-uniform occupation of configu- 
ration space. This alone does not imply non-equilibrium 
behavior in the absence of a DOS (or phase space mea- 
sure) for non-Hamiltonian dynamics [13[. For a beneath 
the critical value a c = |, a region of dynamically inac- 
cessible states appears (see Fig. [I]), which may not be 
explained by a DOS. This is seen if the system is coupled 
to a reservoir, when these states become accessible and 
populated. As a decreases, the inaccessible states take 
over more of configuration space, until at the maximally 
dissipative limit (a = 0) the only occupied states are 
with a single particle holding all the system's energy. 

We now calculate the threshold in a for inaccessible 
states to exist. As a decreases below a c , the first state to 
become inaccessible is that where all particles have equal 



energies. For smaller values of a, the inaccessible region 
grows around this equal energy state. For N = 3 we 
substitute the final energies e[ = e' 2 = e' 3 = (e) = jjf in 
Eq. (pQ), and see that the energy of particle k prior to the 
last interaction had to be = (3 — ~) (e). All energies 
remain positive by the dynamics, hence in order for this 
e/c to be positive, a must be larger than a c = |. We 
therefore conclude that when a < a c the state e\ = e2 = 
es may not be reached by the dynamics, although it is 
energetically allowed. For N > 3 we similarly check what 
is the initial state required for the system to reach e[ = 
. . . = e' N = (e) after a single interaction. Again, Eq. (pQ) 
implies that the energy of particle k before the interaction 
is ek = (3 — ^) (e), leading to the same relation as above. 
To conclude, due to the irreversibility of the dynamics an 
isolated system has inaccessible states for a < a c = ~, 
irrespective of the number of particles in the system. 

To compute the size of the inaccessible hole in configu- 
ration space, we note that for a < a c the accessible state 
closest to the inaccessible equal energy state is that ob- 
tained from the pre-collisional state with e^ = 0, e^+ej = 
3(e), and all other particles having energy (e). The in- 
teraction then brings the system to = ej = ^(e), 
ek = 3(1 — a)(e). This is the nearest the system can get 
to the equal energy state, as measured by the distance: 
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The inset to Fig. [T] verifies this form for the minimal 
distance to the equal energy state for general TV. 

Our dynamics are irreversible, and it is of interest to 
quantify the resulting DB violation, which may serve as a 
measure of how far a system is from equilibrium. To this 
end, we consider a single particle in a large isolated sys- 
tem. Figure [2] depicts the particle's steady state transi- 
tion rates R(e, e')dede' from a region de around e to a re- 
gion de' around e' [14]. In equilibrium these are symmet- 
ric irrespective of the DOS, that is i?(e, e') = R(e' , e) for 
any pair of energies. In our case, i?(e, e') 7^ R(e f , e), and 
we quantify the overall DB violation by the total proba- 




FIG. 2: (Color online) Transition rates R(e,e) of a single 
particle from energy e to energy e for a — 0.4. Colorbar 
indicates the log of the transition rate. 
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FIG. 3: (Color online) a) Total probability current J, tail 
decay constant Q, and normalized second moment Ci vs a. 
b) Single particle energy distribution for various values of the 
restitution coefficient a (indicated next to each line). 



bility current defined as J = J [R(e, e') — R(e! ', e)] 2 dede' 
[Hj]. J vanishes in the equilibrium limit a = 1 and in 
the singular limit a = (where transition rates between 
the extreme states of a single particle possessing all the 
energy are symmetric). It exhibits two peaks with a min- 
imal DB violation around a « 0.6 (see Fig. [3^) [16[. 

Interestingly, the single particle energy distribution 
p(e) exhibits minimal deviation from the Boltzmann dis- 
tribution around a w 0.6 as well. We use the procedure 
described in [13, EH to calculate p(e) in terms of all its 
moments (e n ), and find that each normalized moment 
C n = L may be expressed in terms of lower moments: 
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with K n = ^j + (l — a) n . The distributions for different 
values of a have exponential high energy tails with differ- 
ent decay constants (see Fig. [3)3), but share the same first 
moment (e). As in PJJ], we characterize the distribution's 
high energy tail by ratios of succeeding high order mo- 
ments. For a p(e) ~ exp[— e/(Q(e))] tail, as n — > 00, one 
has C n /C n -i — > Q with Q finite. Iterating Eq. ([3]) we 
find that C n /C n -\ reaches a constant value Q(a) which 
switches from > 1 to < 1 at a « 0.58 (see Fig. [3^). 

A measure for deviation from exponential behavior at 
moderate energies is the normalized second moment 
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C2 switches from > 1 to < 1 at a = 0.6 (see Fig. [3^l) . 
For our model it is equal to the ratio of two effective tem- 
peratures in the system C2 = Tfd/Tq OH EH 03 : the 
fluctuation-dissipation temperature Tfd scales fluctua- 
tions around the steady state to corresponding response 
functions; the granular temperature Tq = (e) measures 
the average energy per degree of freedom. 

We now consider thermal contact between two systems, 
A and B, comprised of particles of different types. We 
assume the irreversible dynamics of Eq. (pQ) with restitu- 
tion coefficients olaa, olbb, an d olab — &ba depending 
on the species of the first two particles (i, j) in each inter- 
action. For example, when an A particle and a B particle 



collide, the restitution coefficient is olab irrespective of 
the type of particle k receiving the dissipated energy. 

We define a dimensionless coupling strength / between 
the systems, so that / = corresponds to no contact 
and / = 1 to intimate contact as in a homogeneous mix- 
ture. We consider the case where the two systems have 
the same number of particles N and exhibit the same 
interaction rate per particle. Therefore, of the eight pos- 
sible ways to choose three particles from the two sys- 
tems, each of the six inter-system interactions (AAB, 
ABA, ABB, BAA, BAB, BBA) occur with probability 
f /8, and each of the two intra-system interactions (AAA, 
BBB) with probability (4 - 3/)/8. We solve the model 
for general / and concentrate here on the weak coupling 
limit < / < 1. 

Averaging the energy transfers between the systems 
over these possible interactions we find that 
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This enables us to identify To = (2 — a)(e) as the oper- 
ational temperature 0, Q of our model: heat flows from 
the system of high To to the one with low To until they 
eventually equalize. Due to the energetic isolation of our 
model, To is defined for the individual systems and, un- 
like open systems 0, depends only on the system's prop- 
erties and not on the contact details (/ and a^). 

To of a non-equilibrium system may be measured by 
attaching to it a thermometer with a = 1. A small 
enough thermometer will not affect the system. After 
the systems reach a mutual steady state we detach the 
thermometer from the system and allow it to reach equi- 
librium. Its (true) temperature then is equal to the oper- 
ational temperature of the measured system since in the 
equilibrium limit a = 1, Tq gives the temperature. 

Other definitions of effective temperatures do not 
equalize upon contact and hence do not behave as ther- 
modynamic temperatures: The granular temperature 
satisfies Tq = (e) — To I (2 — a), and thus differs across 
the systems; The fluctuation-dissipation temperatures in 
each of the two systems in contact may be calculated for 
general coupling, and shown to depend on the coupling 
details. In the weak coupling limit the contact affects 
only the partitioning of energy between the systems and 
not the normalized energy distributions within each of 
them. Therefore, in each system Tfd m &y be evaluated 
from T G by Eq. gj) and using C 2 = T FD /T G . This yields 
different values of Tfd in the two systems. 

In equilibrium, temperature equalization reflects the 
more fundamental entropy maximization. The equaliza- 
tion of an effective temperature To does not necessarily 
imply that entropy is maximized here. Therefore, we 
now test the second law of thermodynamics and consider 
what happens to the total entropy of two systems when 
they are connected. We begin with disconnected systems 
of different types of particles (aAA 7^ &bb), each in its 
own isolated steady state, with some partitioning of en- 
ergy between them. We then connect the two and allow 
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them to reach a mutual steady state. At this time, we 
disconnect them and wait for each one to reach its new 
steady state. Long after they are disconnected, the effect 
of the contact is to have repartitioned the energy without 
modifying the normalized energy distributions in each of 
the systems, and each system's total energy is its only 
energy scale. We may therefore write the multi-particle 
energy distributions in the systems as 

( \ \~ N ( ei eN \ 

p A {e u ...,e N ) (E A ) if I ^"'•••'^"5^AA J 

p B (e u ...,e N ) = (E B )~ N cp (j^,...,^;a BB ^J (6) 

with (f a dimensionless function of the dimensionless en- 
ergies {ei/E} and the dimensionless model parameters 
and a BB in this case). The systems' entropies 
S=— J p(ei, . . . , ejv) lnp(ei, . . . , CN)dc\ • • • dcN now scale 
as Sa = N\u(Ea) + const, S B = N\n(E B ) + const, with 
the additive constants depending only on the normalized 
energy distribution within each system and not on the 
way energy is partitioned between the systems, that is, 
only on a a a and a BB but not on E A and E B jl8j ]. 

When the systems are disconnected entropy is unques- 
tionably additive. Therefore, these expressions for the 



entropy imply that equipartition of energy between the 
systems would maximize the total entropy. Since, how- 
ever, the irreversibility of the dynamics leads to Ea^E b , 
we conclude that entropy is not maximized upon con- 
tact [19]. In particular, if two systems of equal energy 
per degree of freedom are connected, their entropy will 
decrease. That is, unlike the equilibrium case, the sec- 
ond law of thermodynamics does not hold and entropy 
defined in the standard way may decrease rather than 
increase when a constraint is removed [201 ]. 

In conclusion, we characterized non-equilibrium steady 
states of an energy conserving model, and investigated 
what happens when two such systems are brought into 
thermal contact. Although effective temperatures gen- 
erally do not equalize, we identified an operational tem- 
perature which determines heat flow, ultimately reaching 
a common value when mutual steady state is achieved. 
Moreover, we demonstrated explicitly that the total en- 
tropy may decrease in this case. 
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